Do Fermions and Bosons Produce the Same Gravitational Field? 
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We examine some cosmological consequences of gravity coupling with different strength to 
fermions and bosons. We show that this leads to a different perturbation of the standard pic- 
ture of primordial nucleosynthesis than the addition of extra neutrino types or overall scaling of the 
value of G. Observed abundances of deuterium and 4 He place bounds on the ratio of the bosonic 
gravitational constant (Gb) to the fermionic gravitational constant (Gf) of 0.45 < Gb/Gf < 0.92 
at la and 0.33 < Gb/Gf < 1-10 at 2a. A value of Gb < Gf can reconcile the current "ten- 
sion" between the abundances of deuterium and 4 He predicted by primordial nucleosynthesis. We 
comment briefly on other cosmological effects. 



We examine some cosmological consequences of as- 
suming that gravity is not blind to statistics and that 
the gravitational coupling constant, G, is different for 
fermions and bosons. The gravitational constant is noto- 
riously difficult to measure to high precision and is poorly 
known in comparison to other fundamental constants 
Ly. LJ LJ] Hence there has been extensive experimental 
and theoretical exploration of possible deviations from 
the standard Newtonian and general relativistic theory 
of the gravitational interaction |j] . These investigations 
have recently been rejuvenated by the realisation that 
"large" extra dimensions can be accommodated in string 
theories and can lead to anomalous gravitational interac- 
tions on scales as large as 10 -5 m |5|,|fj. Supersymmetry 
assumes that there is a fundamental symmetry between 
fermions and bosons but this symmetry must have been 
broken at the TeV scale. Perhaps this breaking produces 
further asymmetries in fermion-boson properties or cou- 
plings at lower energies? 

In allowing G to be different for bosons and fermions, 
the key question is the scale over which a particle is con- 
sidered a boson or a fermion. We have chosen to set this 
scale equal to the nucleon scale, i.e., all nucleons are con- 
sidered "fermions" in our discussion. There are certainly 
alternatives to this approach. For example, one could 
treat the constituents of the nucleons as the fundamen- 
tal particles, with the quarks coupling to G as fermions 
and the gluons as bosons; this is certainly the correct ap- 
proach prior to the quark-hadron phase transition. Alter- 
natively, one could define bosons and fermions on a larger 
scale, so that a helium nucleus, for instance, would cou- 
ple to gravity as a boson rather than as four fermions. 
We believe that our approach to this admittedly specula- 
tive topic is the most logical, but it is important to note 
that other definitions of bosons and fermions would yield 
other results. 

Note that we assume that this differential coupling of 
bosons and fermions to gravity affects only the source 
term in the Einstein equations. We assume that the 
Equivalence Principle still holds, so that a boson and 
a fermion in a given gravitational field will still follow 



exactly the same trajectories. Because of this distinc- 
tion, the gravitational coupling to bosons is not probed 
experimentally. All direct terrestrial experimental mea- 
surements of G are made with a fermionic source (see ref. 
|4| for a review), so that the measured G is really Gf, 
the coupling to fermions, and the bosonic G, denoted 
Gb, is undetermined. Observations of the gravitational 
deflection of light do not isolate Gb either, since they 
measure the behavior of bosons under a fermionic source, 
and under the assumptions noted above, a boson in the 
gravitational field generated by a fermion will follow the 
same trajectory as a fermion in this gravitational field. 

The Shapiro time-delay effect has been considered as a 
constraint on Equivalence Principle violation under the 
assumption that gravity couples differently to neutrinos 
and photons in the PPN approximation. A comparison 
of the difference in arrival times for neutrinos and pho- 
tons from SN1987A allows a bound to be placed on any 
PPN 7 parameter difference between photons and neu- 
trinos, With \jphoton -lneutrino\ < O(10" 3 ), 0,0, with 

small uncertainties due to the gravitational field of the 
Galaxy. Differences in gravitational coupling to v e and 
v e have also been studied for SN1987A, by Pakvasa et 
al [3, yielding |7„ e — 7p e | < 10~ 6 . It is interesting to 
note that a difference only of order 10 -14 in the coupling 
of gravity to v e and v^ would make the gravitational 
transformation of v e to v^ of similar strength to the con- 
version due to the MSW effect 0.E1E1- In what follows 
we will consider the simplest scenario of a difference in 
gravitational coupling constant for bosons and fermions, 
but assume that the coupling is the same for particles 
and antiparticles and for all fermionic species. These as- 
sumptions can be relaxed straightforwardly if required. 
(See, e.g., Ref. 12] for differential coupling to different 
families, and Ref. [b| for different coupling to particles 
and antiparticles). 

The only place where bosons act as a significant source 
term for gravity is in the overall expansion of the universe 
as a whole, in which photons contribute significantly (at 
least in the early universe) to the overall energy density. 
In seeking to constrain Gb, therefore, the best place to 



look is at early universe tests such as Big Bang Nucle- 
osynthesis (BBN) and the cosmic microwave background 
(CMB). We will explore the former in detail and com- 
ment briefly on the latter. 

In the standard flat Friedmann cosmological model, 
the overall Hubble expansion rate is given at early times 
by 



H 



-nGp 



1/2 



(1) 



where p is the total gravitating density. Now consider a 
model in which we have two separate gravitational con- 
stants: Gf for fermions and Gb for bosons. In our model, 
Gf is just set equal to the measured value of Newton's 
constant, G = 6.6742 ± 0.0010 x 10~ 8 cm 3 gm" 1 s" 2 Q, 
so all of the information in the model is contained in the 
ratio Gb/Gf, which we will denote by /bf'- 



Jbf 



Gl 
G, 



Then equation Q for the expansion rate becomes 
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H ~ ( -tt(GbPb + GfPf] 
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and the problem is equivalent to changing the density 
of bosons by the factor of fsF and leaving the fermion 
density unchanged. 

This model is qualitatively similar to both adding (or 
subtracting) relativistic degrees of freedom, or to chang- 
ing the overall value of G, both of which have been exten- 
sively explored in connection with BBN, but we will now 
show that it is different from either of these previously 
studied situations. 

Consider first the effect on the expansion rate of adding 
additional relativistic degrees of freedom to p. (The 
resulting variation in the primordial element produc- 
tion can be used to constrain such a change, e.g., Refs. 
[H, EH EH H3Hi Ull)- It is conventional to parametrize 
such a change in terms of the effective number of addi- 
tional (i.e., beyond v e ^v ili ^v T ) relativistic two-component 
neutrino degrees of freedom, AN V . In this case the total 
relativistic energy density prior to e + e~ annihilation is 



P = [2 + (7/8)10 + (7/8)2AAy _ T 4 



(5) 



where the first term in square brackets counts the boson 
degrees of the freedom (photons), the second counts the 
fcrmionic degrees of freedom (e + e~ and vv), the third 
counts any hypothetical additional relativistic degrees of 
freedom, and T is the photon temperature. (We use units 
with H = c = kg = 1 throughout). After e + e~ annihila- 
tion, when the photons are heated relative to the neutri- 



nos, the corresponding energy density is 

P = [2 + (7/8)(4/ll) 4 / 3 6 + (7/8)(4/ll) 4 / 3 2A^]^-T 4 . 

These expressions can be rewritten in terms of the photon 
density, p 7 , to give [l8| 

Pbefore = 5.375[1 + 0.1628 AAyp 7 , (7) 

p af ter = 1.681 [1 + 0.1351 AJV„]p 7 , (8) 

where the subscripts before and after refer to the den- 
sity before and after e + e~ annihilation, respectively. 

Another way to parametrize the change in the expan- 
sion rate is through a change in the overall value of G, or, 
equivalently, multiplying equation P by a "speed-up" 
factor S. By incorporating such a change into BBN, the 
resulting changes in the element abundances can be used 
to constrain a time-shift in the value of G J20l l2ll l22l |23J| . 
As emphasized by Kneller and Steigman (18| . these two 
ways of modifying BBN (adding additional relativistic 
degrees of freedom, or multiplying G by a constant) while 
qualitatively similar, are inequivalent. This can be most 
easily seen by noting that changing G by some fixed fac- 
tor fc is completely equivalent to changing p by this 
same factor; the "effective" p which enters into equation 
(Q is then just 



Pbefore = 5.375 [/g]P 7 , 
Pafter = l.68l[fc]P~ r 



(9) 
(10) 



In order to make a change in G equivalent to a change in 
the number of effective neutrino degrees of freedom, we 
would need fa = 1 + 0.1628AA„ before e + e~ annihila- 
tion and fa = 1 + 0.1351AA^ after e + e~ annihilation; 
obviously, it is impossible to satisfy both equations si- 
multaneously. 

Now consider the equivalent expressions in our model, 
when Gb ^ Gf- In equations (|5I6|I . the first term in 
brackets gives the contribution to the energy density from 
the photons, which are the only bosonic degrees of free- 
dom present during BBN. Hence, when Gb ^ Gf, the 
effective density becomes 



Pbefore = [2/ BF + (7/8) 10] — T 4 , 



(11) 



Pafter = [VbF + (7/8) (4/1 1) 4 / 3 6] ^T 4 . (12) 

These expressions can be rewritten in terms of the photon 
density as 



Pbefore = [fBF + 4.375]p 7 , 
Pafter = UbF + 0.681]p 7 . 



(13) 
(14) 



A comparison of equations (|13fl - (|14l) with equations (J7J) 
- © and with equations 10 - (|10fl shows that changing 
Gb/Gf is inequivalent (in terms of its effect on the ex- 
pansion rate) to changing either the overall value of G, 



or adding additional relativistic degrees of freedom. This 
can be seen most easily by fixing pbefore to the same 
multiple of p 7 in all three cases; it is easy to see that the 
resulting values of p a fter are all different. 

We now consider the effect of taking Jbf 7^ 1 on the 
primordial element abundances. The primordial produc- 
tion of 4 He is controlled by the competition between the 
expansion rate and the rates for the weak interactions 
which govern the interconversion of neutrons and pro- 
tons: 



Using a modified version of the Kawano nucleosynthe- 
sis code J2g we scan over the (77, Jbf) plane, calculating 
the likelihood for a given pair of values. (In doing so, 
we take the distribution of \og(D/H) to be Gaussian as 
in, e.g., Refs. |23|,|25j, rather than taking the distribu- 
tion of D/H to be Gaussian as in, e.g., Ref. |l]|. In 
practice, this should have only a small effect). The 1 — a 
and 2 — a contours are shown in Fig. 1 (where we take 
i] = 7?io x 10~ 10 ). The limits on Jbf, at the 1 — a and 



n + v e «-> p + e , 
11 + e + «-» p + v ei 

n <-> p + e~ + v e . 



(15) 



At high temperatures, T > 1 MeV, the weak-interaction 
rates are faster than the expansion rate, H, and the 
neutron-to-proton ratio (n/p) tracks its equilibrium value 
exp[— Am/T], where Am is the neutron-proton mass dif- 
ference. As the universe expands and cools, the expan- 
sion rate becomes too fast for the kinetic equilibrium to 
be maintained by weak interactions and n/p freezes out. 
Nearly all the neutrons which survive this freeze-out are 
converted into 4 He as soon as deuterium becomes stable 
against photodisintegration, but trace amounts of other 
elements are produced, including deuterium (see, e.g., 
Ref. [24( for a review). Therefore, the primordial pro- 
duction of 4 He is very sensitive to the expansion rate of 
the Universe at temperatures ~1 MeV, so BBN has been 
used many times to constrain any change in this expan- 
sion rate. 

As is the case for other models which change the expan- 
sion rate, the primordial deuterium abundance is most 
sensitive to changing the baryon-photon ratio 77, and it 
essentially provides the upper and lower bounds on 7;. 
The predicted abundance of 4 He within this range for 77 
can then be calculated as Jbf varies, allowing bounds to 
be placed on Jbf- (This is something of an oversimplifi- 
cation, as the deuterium abundance also depends weakly 
on Jbf] our calculation correctly incorporates this de- 
pendence) . 

The primordial abundance of deuterium has been in- 
ferred from QSO absorption systems. We use the abun- 
dance estimated by Kirkman et al. 25]: 



log{D/H) = -4.556 ± 0.064, 



(16) 



where all errors are quoted at the 1 — a level. 

The abundance of 4 He can be inferred from low- 
metallicity HII regions, but there are significant discrep- 
ancies between different estimates (see Ref.pjj for a re- 
cent analysis, and references therein). We will follow Ref. 
[ijj and take the primordial 4 He mass fraction, Yp, to be 



Y P = 0.238 ±0.005, 



(17) 



an estimate which is consistent, for example, with that 
found in the review of Olive et al. [2J| . 
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FIG. 1: The 1 — a (solid) and 2 — a (dashed) contours in the 
7710, Jbf plane, where 77 = 7710 x 10 -10 and Jbf — Gb/Gf- 



2 — a levels, are: 



0.45 < f BF < 0.92, (1 - a), 
0.33 < Jbf < 1.10, (2 - a). 



(18) 
(19) 



These limits are the main result of our paper. 

Note that the standard model is excluded at the 1 — a 
level. This is not surprising; it is due to the tension 
which currently exists between the observed deuterium 
abundance and the observed 4 He abundance; the former 
prefers a higher value of 77 than does the latter. This sort 
of result, then, also shows up in discussions of changing 
AN V [ijj and in analyses of other allowed deviations from 
standard physics (e.g., Ref. 27]). Note that a value of 
Gb less than Gf provides yet another mechanism for 
resolving this tension. Clearly, Gb 3> Gf is excluded by 
BBN, but Gb could be much smaller than Gf- 

Of course, these conclusions are strongly dependent 
on the parameters that go into the calculation. A higher 
primordial 4 Hc abundance, such as that claimed, for in- 
stance, in Ref. 28], would eliminate the tension between 
the deuterium and He abundances, and shift our allowed 
region upward in Fig. 1. A similar effect would occur if 
the neutron lifetime were measured to be shorter than 
the currently accepted value 29]. 



Including CMB data in our analysis will improve these 
limits, but not by much. Although, as we have empha- 
sized, our model differs both from changing AiV„, or from 
changing G, it is qualitatively similar to both of these. 
For the case of changing AiV„ , the addition of CMB data 
constrains 77 more tightly than BBN alone, but it has 
only a small effect on the overall limits on AiV„, com- 
pared to using BBN alone |l9j. Similarly, BBN provides 
a stronger constraint on a change in G than does the 
CMB (compare, for example, the results in Ref. Q for 
BBN with those in Ref. $jfr for the CMB). 

It is difficult to imagine other environments in which 
the value of Gb could manifest itself. If a scenario of 
baryon number generation by out of equilibrium decay 
of superheavy vector or scalar bosons [3lL l32j | were pre- 
cisely established as the source of the value of 77, then 
there would be a dependence of r\ on Jbf via the ra- 
tio of boson decay rate to the total expansion rate, but 
this is not the case at present. If the dark matter were 
bosonic, then the overall expansion rate would also be 
altered in the matter-dominated era. This could only 



be detected, however, if the number density and mass of 
the dark-matter particle were independently determined, 
rather than being inferred from the Hubble expansion it- 
self. The effect on large-scale structure in this case is 
equivalent to a model with one gravitational interaction 
for dark matter and a different one for luminous matter 
\,Y.'\ . It remains to be seen whether there are further con- 
sequences of Gb 7^ Gf at very early times (when T > > 1 
MeV) which could lead to observable consequences at low 
energy. 
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